A self-consistent theory of atomic Fermi gases with a Feshbach resonance at the 

superfluid transition 
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A self-consistent theory is derived to describe the BCS-BEC crossover for a strongly interacting 
Fermi gas with a Feshbach resonance. In the theory the fluctuation of the dressed molecules, 
consisting of both preformed Cooper-pairs and "bare" Feshbach molecules, has been included within 
a self-consistent T-matrix approximation, beyond the Nozieres and Schmitt-Rink strategy considered 
by Ohashi and Griffin. The resulting self-consistent equations are solved numerically to investigate 
the normal state properties of the crossover at various resonance widths. It is found that the 
superfluid transition temperature T c increases monotonically at all widths as the effective interaction 
between atoms becomes more attractive. Furthermore, a residue factor Z m of the molecule's Green 
function and a complex effective mass have been determined, to characterize the fraction and lifetime 
of Feshbach molecules at T c . Our many-body calculations of Z m agree qualitatively well with the 
recent measurments on the gas of 6 Li atoms near the broad resonance at 834 Gauss. The crossover 
from narrow to broad resonances has also been studied. 

PACS numbers: 03.75.Ss, 05.30.Fk, 74.20.-z 



Introduction. — The recently demonstrated strongly 
interacting Fermi gases of 40 K and 6 Li near a Feshbach 
resonance have attracted a lot of attention 0, 0, IE IE 
S 0| In these systems, because of the coupling be- 
tween a pair of atoms in an open channel and a bound 
molecular state of a closed channel, it is possible to 
obtain an arbitrarily strong attractions between atoms 
by tuning the energy of the molecular state across the 
resonance with a magnetic field. It therefore provides 
us a unique opportunity to experimentally examine the 
crossover from a Bardeen-Cooper-Schrieffer (BCS) super- 
fluidity to a Bose-Einstein condensate (BEC). 

The use of the Feshbach resonance, however, signif- 
icantly complicates the theoretical analysis of current 
experimental results, due to the presence of the bound 
molecular state. It has been argued that a two-channel 
(resonance) model that describes explicitly both the open 
channel and closed channel physics should be adopted as 
a minimum ingredient for the theory [E ^J, 0] . In the 
model, the molecular bound state in the closed channel is 
treated as a featureless bosonic particle, and two essential 
parameters, the threshold energy of molecules and the 
atom- molecule coupling constant, are required to charac- 
terize the total system. This is in contrast to the tradi- 
tional single-channel model, in which the closed channel 
or the "bare" Feshbach molecules are neglected. Thereby 
the model is characterized solely by a single parameter 
of s-wave scattering length a, i.e., the interatomic inter- 
action is described by a contact potential with strength 
A-Kh 2 a/m. The applicability of these models to current 
experiments, together with their equivalence in certain 
limits, are actively debated at present. A closely related 
subtle problem is the nature of Fermi superfluidity near 
the resonance, i.e., whether the superfluid transition is 
associated with small pairs that are molecular in charac- 
ter, or large objects akin to Cooper pairs. 



In this paper we present a many-body theory of the 
two-channel model for an interacting ultra-cold Fermi gas 
with a Feshbach resonance at and above the superfluid 
transition temperature T c . In the theory the strong fluc- 
tuations of the preformed Cooper-pairs and of the Fes- 
hbach molecules have been incorporated within a self- 
consistent T-matrix approximation, similar to that con- 
sidered by Haussmann for single-channel model 0|. We 
then apply the theory to explore the normal phase of the 
gas at BCS-BEC crossover, emphasizing the importance 
of the resonance width. The purpose of the exploration 
is three-fold. (1) First, a reliable T c is predicted. Till 
now there are strong evidences for the emergence of su- 
perfluidity at the crossover 0, [E IE 0, IE IE| • Precise ex- 
perimental characterizations of the thermodynamic and 
dynamical properties of strongly interacting Fermi gases 
are also undergoing. In this respect, a reliable theoretical 
prediction of T c of the crossover gas is, therefore, useful 
for guiding future experimental investigations. Previous 
calculations of T c is based on the so-called Nozieres and 
Schmitt-Rink (NSR) scheme EE 03, which is an anal- 
ysis equivalent to the ladder approximation using free 
fermionic Green's functions. In both weak and strong 
coupling limits this approximation leads to correct re- 
sults [IjI . However, at the crossover it becomes less ac- 
curate as the fermionic quasiparticles are far from being 
free particles. Our theory has the advantage of taking a 
self-consistently dressed Green function in the T-matrix, 
and is expected to give a reasonable T c than the NSR 
scheme. (2) Second, to characterize the nature of pairs 
at the crossover, the fraction and lifetime of the "bare" 
Feshbach molecules have been studied, by extracting the 
residue factor Z m and the effective mass from the Green 
function of molecules. The value of Z m can be directly 
determined in experiments using the Stern-Gerlach selec- 
tion technique to measure the magnetic moment of the 



pairs [15], or using the optical spectroscopy to project 
out the singlet component of the pairs [16]. Our cal- 
culations of Z m emphasize the many-body effects at the 
crossover, in particular, on the BCS side of the resonance. 
For an ultracold gas of 6 Li atoms at the broad resonance 
B = 834 G, we find a reasonable agreement between our 
theoretical predictions and the recent spectroscopy mea- 
surements by Hulet et al. (3) Finally, the crossover 
from a narrow to a broad resonance with increasing atom- 
molecule coupling strengths has been investigated. 

Self-consistent T-matrix approximation. — We con- 
sider a homogeneous gas of fermionic atoms in two dif- 
ferent hyperfine states in the vicinity of a Feshbach res- 
onance. In the grand canonical ensemble the gas of 
N atoms is described by an atom-molecule Hamiltonian 

H = ( £k ~ V) C £x C k<x + U bare Y C k+qT C k'-ql C k'lC kT 

kcr kk'q 

q 

+9bare ^ ( fo q c k+q/2T c -k+q/21 + k.C.) . (1) 
kq 

Here the two hyperfine states are denoted as pseudospins 
a =\,\, with N-\ — N\ — N/2. c^ a and 6+ represent 
the creation operators of an atom with the kinetic en- 
ergy ek = ft 2 k 2 /2m and of a Feshbach molecule with the 
dispersion e q /2 + v ba re = fi, 2 q 2 /4m + Vbare, respectively. 
The magnetic detuning or the threshold energy Vbare is 
a key parameter in this model, and can be conveniently 
adjusted by an external magnetic field B. It changes sign 
across the resonance from above and enables the forma- 
tion of stable molecules. As the molecule is made of two 
atoms, a chemical potential \i is introduced to impose the 
conservation of the total number of the bare Fermi atoms, 

i.e., (AO = (Ek. 4a c ^) + (Eq&J&q) = N - Moreover, 
two contact interactions have been adopted for the non- 
resonant (or back-ground) atom-atom scatterings and the 
atom-molecule coupling, as parameterized by Ubare and 
gbare, respectively. 

As usual, the use of contact interactions leads to an 
(unphysical) ultraviolet divergence, and requires a renor- 
malization that expresses the bare parameters Ubare, 
gbare, and i^bare, in terms of the observed or renormal- 
ized values Uq = A-n:h 2 abg/m, go, and Vq = A/i (B — Bq), 
where a bg is the back-ground s-wave scattering length for 
the atoms, go is associated with the width of resonance, 
and A/j, is the magnetic moment difference between the 
atomic and the bound molecular states. The renormal- 
ization can be done by solving either the Lippmann- 
Schwinger scattering equations or the vertex functions 
in the two-body limit. The outcome is very transpar- 
ent in physics: one can just replace in everywhere the 
bare parameters by their corresponding observable val- 
ues. In addition, the divergent Cooper pair or two- 
particle propagator (defined later) could be regularized 
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Figure 1: Diagrammatic representations of the atomic self- 
energy (a), the molecular self-energy (b), the Bethe-Salpeter 
equation for the four-point vertex function (c), and the atom- 
molecule vertex function (d). The solid (dot-dashed) line rep- 
resents G (Do), the dotted line describe Ubare, and the empty 
circle is the bare atom-molecule coupling gbare- 



as n(q, iv n ) -> n(q, %v n ) - £ k 1/ (2e k ). In Ref. [H, it 
is stated that this procedure yields the correct expression 
for the molecular binding energy, and therefore incorpo- 
rates the exact two-body scattering process. A similar 
conclusion has also been arrived by Falco and Stoof (lflj . 

To relate the value of go to experimentally known 
parameters, we compare the expression for the ef- 
fective interactions between atoms in the low-density 
limit, U e ff = Uo + go/ (~^o) 5 with the s-wave scat- 
tering length of the atoms as a function of magnetic 
field, a(B) = a bg [1 + AB/(B — B )], where AB is 
the width of resonance. Using the relation U e ff = 
ATrh 2 a(B)/m and eliminating Uq and vq, we have that 
go = 2fi [n (—abgABAfi) /m] 1 ^ 2 . For the 6 Li atoms at the 
broad resonance Bo — 834 G of interest here, AB 300 
G, Afi w 2fiB, and a bg ~ — 1405ao 0], where \xb is the 
Bohr magneton and ao — 0.529 A is the Bohr radius. 

We solve the atom-molecule Hamiltonian in the 
self-consistent T -matrix approximation by summing up 
infinite ladder diagrams involving Ubare and particle- 
particle bubble diagrams involving gbare- We consider 
here only the normal phase of the atomic gas. The gen- 
eralization to the superfluid phase is more interesting 
(and, of course, more complicated), and will be addressed 
elsewhere. As shown in Fig. la, diagrammatically the 
dressed Green function of atoms G(k, iu m ) and its self- 
energy £ a (k, iui m ) read, respectively, 

G~ x (k.,iu m ) = Gq 1 (k, iuj m ) - S (k, iu m ), (2) 
S a (k,io; m ) = T T(q, iv n )G (q - k, iv n - iu m \$) 



where Go (k, iu> m ) = [iu> m — (e^ — fi)] 1 is free atomic 
Green function, and ius m (iv n ) is the fermionic (bosonic) 
Matsubara frequencies. The approximate four-point ver- 
tex function r(q, iv n ) is solved readily from the Bethe- 
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Salpeter equation (see Fig. lc), 

r _1 (q, iv n ) = U~^ f (q, iv n ) + II(q, iv n ), (4) 

where II = TJ2k^ m G (q - k, iv n — iui m ) G (k, iu m ) — 
J2k 1/ (2ek), and the effective pairing interaction U e ff = 
Uq + g 2 D a (q, ivn) consists of the background inter- 
action U{) and the induced interaction mediated by 
Feshbach molecules %Dq (q, iv n ), where Do(q,iv n ) — 
1/ [ivn — £ q/2 — + 2/i] is the non-interacting molecu- 
lar Green function. Note that the induced interaction 
dominates near the Feshbach resonance, and acquires a 
dependence on the relative energy of the colliding atoms 
(i.e., retarded effects) that becomes stronger for a nar- 
row resonance. Analogously, it is straightforward to write 
down the atom-molecule vertex function, 



9 1 (q, Wn 



90 1 



9q 1 U U(q,iu n ), 



(5) 



whose diagrammatic equivalent is shown in Fig. Id. This 
leads to the Dyson equation for molecules (see Fig. lb), 

D~ x (q,w n ) = D 1 (q, iv n ) - S m (q, iu n ), (6) 
£ TO (q,iz/ n ) = -goIL(q, iv n )g(q, iu n ), (7) 

where -D (q, iv n ) and S m (q,w„) are, respectively, the 
molecular Green function and self-energy. 

Eqs. 0-0 form a close set of equations, which we 
refer to as a self-consistent theory for the two-channel 
model. We emphasize that in these equations, we have 
used a fully dressed fermionic Green function that has to 
be determined self-consistently. This is contrasted to the 
NSR scheme considered by Ohashi and Griffin, where 
a free fermionic Green function is adopted [1J|. Actu- 
ally, with replacing G by Go in II and S a and writing 
G = Gq + GoS a Go, we reproduce the NSR formalism. 
Furthermore, if we set U e ff = Uq in the four-point vertex 
function and neglect the molecular part, we recover 
the self-consistent theory for the sing le-channel model, as 
proposed earlier by Haussmann fl2|. 

We have solved Eqs. 0-0 in a self-consistent manner 
by an iteration procedure, which is repeated successively 
to obtain G and D in higher orders until convergence is 
achieved. At a finite temperature T > T c , the only quan- 
tity to be determined is the chemical potential /i, which 
has to be adjusted to satisfy the constraint for the total 
density n to t — nF + 2n m , where np — 2G(x = 0, r = _ ) 
and n m = —D(x = 0, r = 0~) are the density of atoms 
and molecules, respectively. In the calculation, the en- 
ergy and length are taken in units of €f = h 2 k F /2m 

and hp 1 , where kp — (37r 2 n tot ) 1 ^ 3 is the characteristic 
Fermi wavelength. Without further clarity, in the calcu- 
lation we take a weak back-ground atom-atom interac- 
tion, kpa bn = —0.1, in accord with the real experimental 
situation jig . 

It is worth noticing that the properties of Fes- 
hbach resonaces is naturally characterized by a di- 
mensionless atom-molecule coupling constant g = 
goh~ 2 m (2-KkF ) X ^ 2 , quoted later as an intrinsic width 
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Figure 2: (color online) (a) The superfluid transition tem- 
perature T c as a function of the effective coupling strength 
v e ff for dimensionless atom-molecule couplings g = 10, 3, 
and 1. The full lines with solid symbols represent the numeri- 
cal results of our self-consistent theory. The dotted lines with 
empty symbols are obtained by NSR method, (b) The popula- 
tion of bare Feshbach molecules n m at T c , and the population 
of preformed Cooper pairs 2n cp = uf — «fo at T c , where 
Ufo = 2 / (ek — m), and f(x) is the Fermi distribution 
function. 



of resonances. To illustrate the physical implica- 
tion of this dimensionless parameter, we note that 
the important energy scale of resonances can be de- 
fined as e R = A/it (AB) m , where (AB) m /AB = 



(AfiAB) I [h 2 / ( 



[2il l2j- g is related to the 
ratio between en and the Fermi energy eF via g — 
[eji/ep] 1 ^ 4 ■ Hence, if the Fermi energy is well within 
the characteristic energy scale of resonances, i.e., g 1, 
we have a broad resonance. The opposite limit of g < 1 
defines a narrow one. We note also that it is possible 
to tune the value of the intrinsic width of resonances by 
changing the density of the gas, as g oc kp 1 ^ 2 oc n tc ^/ 6 . 

Superfluid transition temperature T c . — The super- 
fluid transition temperature can be conveniently deter- 
mined by the Thouless criterion, 



[r(q = 0,iz/ n = 0)] T * =0, 



(8) 



which describes the instability of the normal phase of the 
gas towards the formation of Cooper pairs of atoms [lflj . 
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In our T-matrix approximation, the Thouless criterion is 
consistent with the Hugenholtz-Pines condition for Fes- 
hbach molecules that the molecular Green-function (jjjj) 
must develop a pole at q =iv n = as T — > T c . As the 
four-point vertex function represents roughly the propa- 
gator of Cooper pairs, this consistency implies that the 
superfluid phase transition is accomplished by dressed 
molecules that involve both the preformed Cooper pairs 
and the Feshbach molecules (T3 |. 

In Fig. 2a we present our results for T c as a function 
of an effective coupling constant v e ff = 1/ (kp a e //) as 
full symbols, where 47rS 2 a e ///m = Uq + g^/ (— vq + 2fi). 
For comparison the NSR predictions are also plotted by 
empty symbols. Three values of dimensionless atom- 
molecule coupling g have been considered, corresponding 
to the situations with broad (g = 10), medium (g = 3), 
and narrow (g — 1) resonances. For all these intrin- 
sic resonance widths, it is obvious that our results of T c 
increase monotonically with increasing the effective cou- 
pling constant. In contrast, for the broad or medium res- 
onance, the NSR theory predicts a slight hump structure 
at the crossover. This maximum presumably is an arti- 
fact of the lack of self-consistency. As stated by Hauss- 
mann [12}, a repulsive interaction between the pairs leads 
to a depression of T c due to the composite nature of 
pairs. The approximating use of the free fermionic Green 
function in NSR scheme neglects the possible repulsion 
between dressed molecules, and therefore overestimates 
the transition temperature. Nevertheless, in two limits 
of weak and strong interactions, both theories give al- 
most the same results: In the weak-coupling limit of 
v e ff —>■ — oo, they give back the standard BCS result 
of ksTBcs/^F — > 0.614 exp [7rw e ///2], as shown by the 
dashed line; On the other hand, in the strong-coupling 
limit of v e ff — > +00, they correctly reproduce the BEC 
temperature of kgTgEC / e F ~ 0.218. 

The transition temperature depends also on the in- 
trinsic width of Feshbach resonances in a non-trivial 
way. As shown in Fig. 2b, theoretically, with decreas- 
ing the intrinsic width the character of dressed molecules 
changes from the large object of Cooper pairs to the 
small point structure of "bare" Feshbach molecules. As 
a fact that the transition temperature of non-interacting 
"bare" molecules is larger than that of performed Cooper 
pairs, the less population of Cooper pairs tends to en- 
hance our predicted transition temperature. This feature 
should be universal with respect to the inclusion of repul- 
sive molecule-molecule interactions in the atom-molecule 
model Ql, because, contrasted to the composite Cooper 
pairs, the transition temperature of a weakly interacting 
gas of true molecules increases with increasing repulsions 
between moleclues jl^l , and thereby is always larger than 
that of Cooper pairs. We note that, interestingly, the 
NSR results go exactly in the opposite direction as the 
intrinsic width decreases, which is again due to the lack 
of self-consistency. 

Molecular probe of the BCS-BEC crossover. — At 
the crossover the dressed molecules are extremely weakly 
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Figure 3: (color online) (a) The real part of Z m at the su- 
perfluid transition as a function of v e ff. (b) ReZ m (T c ) and 
Re Z m (2T C ) at the crossover for g = 10 in logarithmic scale, 
(c) Temperature dependence of ReZ m at w e // = — 1, 0, and 
+1 for g = 10 . 



bound objects, depending strongly on many-body ef- 
fects. They can hardly be detected by the standard 
technique such as the time-of-flight expansion. In con- 
trast, the "bare" Feshbach molecules are well defined 
point structures, and can easily be characterized by the 
usual optical spectroscopy methods. Therefore it is pos- 
sible to probe the crossover indirectly through the mea- 
surement of "bare" molecules. As the dressed molecules 
contain both Cooper pairs and "bare" molecules 0], i.e., 
| dressed) = e^^/l — Z m |open) + \JZ m | closed), a candi- 
date observable would be the amplitude Z m of the bare 
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Figure 4: (color online) Comparison between our theoretical 
predictions and the experimental data (Fig. 5 in Ref. 
for the residue of the pole of the 6 Li Feshbach molecules at 
the broad resonance Bo = 834 G. 



molecular state in the closed channel. For an ultracold 
gas of 6 Li atoms, the closed (open) channel, to an ex- 
cellent approximation, corresponds to the atomic singlet 
(triplet) state. As a result, the value of Z m can directly 
be determined by selectively projecting out the singlet 
component of the dressed molecules with optical spec- 
troscopy. This is exactly what was done recently by Hulet 
et al. [16], where a tiny value of Z m (~ 10~ 4 ) has been 
observed in the unitary limit (kFCteff — > oo) right on 
resonance. 

In our theory Z m is associated with the residue of the 
pole of the molecular Green function (JJ, which, in the 
spirit of Landau quasiparticle concept, takes the form, 



D (q,ii/ n ) = 



iv n - (h 2 a^/2m* m - fi m ) ' 



(9) 



for sufficiently small q and iv n . Here m* m is a complex 
effective mass, and \i m is a molecular chemical potential 
that goes to zero as T — > T c . In Fig. 3a, we report 
our self-consistent calculations for ReZ m at T c against 
the coupling v e .tu together with the two-body results 
in dashed lines (l7l llS| . In the strong coupling BEC 
region {v e ff > 1) these results are identical as expected. 
However, near the Feshbach resonance they do show a 
large difference, which should attribute to the many-body 
effects. 

The value of Z m decreases significantly with increas- 
ing the atom-molecule coupling. In large g limit where 
the Cooper pairs dominate at the crossover, we find ap- 
proximately Z m (g,v e ff,T) w Z m (v ef f,T)/g 2 . For 6 Li 
gas at the broad resonance Bo = 834 G, we estimate 
g w 120 [ifj. Therefore, from our data at g = 10, the 
expected value of Z m in the unitary limit is in the order 
of 10~ 4 , which agrees qualitatively well with the exper- 
imental findings by Hulet et al. [lfl |. A full comparison 
between our theorectial results and experimental data 
over the entire crossover regime is shown in Fig. 4. We 
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Figure 5: (color online) The real part (a) and imaginary part 
(b) of the effective mass of the bare molecules at the crossover. 



note that this qualitative agreement is robust against the 
unknown temperature in experiment. In Figs. 3b and 3c, 
we show the finite temperature effect on Z m in a logarith- 
mic scale for g = 10. As the temperature doubles, Z m 
increases, but its order is essentially unchanged. We note 
also that in our calculations we do not consider the effect 
of harmonic trap present in experiments, which tends to 
quantitatively decrease the value of Z m . 

At the crossover, the Feshbach molecules can transfer 
or break up into free fermions. This mechanics causes 
the finite lifetime of the molecules, as evidenced by the 
nonzero imaginary part of the residue Z m , or of the effec- 
tive mass m* m . We present in Fig. 5 the results for m* m 
at T c . In the deep BEC limit, the molecules are stable 
so that the imaginary part of m* m is nearly zero and the 
real part Rem*„ « 2m. While in the opposite limit of 
a weak coupling, i.e., v e ff < —1, the imaginary part is 
much larger than the real part, which suggests that the 
molecules are overdamped. In addition, as the dimen- 
sionless atom- molecule coupling g decreases, the value of 
Imm^j becomes smaller, as expected. 

Crossover from the narrow to broad resonances. — 
It was suggested that for a broad Feshbach resonance, 
the gas should be well described by the single-channel 
model with a single parameter (i.e., the s-wave scatter- 
ing length) [23|. This suggestion was supported from 
a coupled-channel calculation in the context of two-body 
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Figure 6: (color online) Illustration of the crossover from nar- 
row to broad resonances on resonance vo = 0. (a) The transi- 
tion temperature T c , (b) the chemical potential at T c , (c) the 
fractions of particles, and (d) Re Z m (T c ) as a function of the 
dimensionless atom-molecule coupling g. 

physics 1231 . or. from a mean-field calculation at zero tem- 
perature j2j. However, the situation is less clear beyond 
the mean-field approximation. In Fig. 6, we show our 
many-body results on resonance (uq = 0) for the critical 
temperature T c , the chemical potential at T c , the residue 
of the pole of "bare" Feshbach molecules Z m , as well as 



the fractions of particles. As g increases, T c and \i (T c ) 
come to saturate, and Z m and n m drop as l/g 2 . There- 
fore, using our self-consist theory we verified numerically 
that in the broad resonance limit the behavior of the gas 
indeed becomes universal, towards the prediction of the 
single-channel model. 

A possible candidate to observe the crossover from nar- 
row to broad resonances is a gas of 40 K atoms at Bo = 202 
G, where we find g ~ 10 under current experimental sit- 
uation jjj. According to the relation g oc n -1 / 6 , it is 
possible to obtain a narrow resonance by increasing the 
center density of the gas, n. This might be achieved by 
using more tight harmonic traps or trapping more atoms 
in the experiment. 

Conclusions. — In this paper we have developed a 
many-body theory of the two-channel model for an inter- 
acting Fermi gas with a Feshbach resonance. The theory 
includes the strong fluctuations in a self-consistent man- 
ner that is necessary at the BCS-BEC crossover. As an 
application, we have studied in detail the normal phase 
properties of the crossover, with a special attention on 
the effects of intrinsic resonance widths. A reasonable su- 
perfluid transition temperature has been predicted. Our 
predictions for the residue of the pole of "bare" molecules 
agree qualitatively well with the experimental findings. 
We expect, therefore, that the present theory gives a re- 
liable description of current BCS-BEC crossover physics 
in atomic normal Fermi gases. The generalization of our 
theory to the superfluid phase is under investigation, and 
will be reported elsewhere. 

We are indebted to Prof. Peter D. Drummond and 
Prof. Allan Griffin for useful discussions. We gratefully 
acknowledge the Australian Research Council for the sup- 
port of this work. 
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